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Abstract 


This  report  summarizes  progress  achieved  this  year  both  in  the  more  general 
formulation  of  our  new  criterion  for  nonlinear  systems  which  allows  us  to  tell  right  away  if  a 
chaotic  system  will  exhibit  a  1/f  spectrum,  and  in  the  application  and  further  study  of  the 
quantum  1/f  effect.  The  general  criterion  was  applied  to  a  one-dimensional  crystal  with 
anharmonic  interactions,  predicting  a  1/f  phonon  number  spectrum  in  the  chaotic  regime  at 
very  low  frequencies  and  always  when  cubic  terms  are  dominant  in  the  potential  energy.  The 
quantum  1/f  theory  was  directly  applied  to  a  quartz  resonator  for  the  first  time,  providing  an 
explanation  for  the  observed  quantum  1/f  frequency  fluctuations  and  optimization  means.  Our 
new  formula  for  collector  1/f  noise  in  ultrasmall  BJT’s  was  found  to  agree  reasonably  with  the 
experiment.  Finally,  the  fractional  dimension  of  band-limited  quantum  1/f  noise  was 
determined  for  the  first  time  and  found  to  coincide  with  the  number  of  octaves  considered,  in 
agreement  with  known  experiments.  This  indicates  a  quantum  chaos  nature  of  1/f  noise  in 
infrared  detectors. 
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I.  INTRODUCTION 


Progress  has  been  achieved  this  year  in  the  study  of  nonlinear  systems  which 
generate  chaotic  1/f  fluctuations,  in  the  application  of  the  Quantum  1/f  Theory  to 
various  materials  used  in  small  and  ultrasmall  electronic  devices,  and  in  the  application 
of  the  Quantum  1/f  Theory  to  electronic  devices.  During  this  year  four  main 
achievements  are  reported.  They  are  our  discovery  of  the  general  mathematical 
principle  causing  the  ubiquitous  1/f  spectrum  in  nature,  science  and  technology  briefly 
presented  in  Sec.  II,  our  application  of  this  principle  as  a  sufficient  1/f  criterion  to 
various  chaotic  nonlinear  systems  such  as  the  one-dimensional  crystal  with  anharmonic 
interactions  presented  as  an  example  in  Sec.  Ill,  our  development  of  a  new  method 
allowing  for  direct  application  of  the  quantum  1/f  principle  to  quartz  resonators  and  to 
other  piezoelectric  and  ferroelectric  systems  presented  in  Sec.  IV,  a  comparison  of  our 
improvements  in  the  quantum  1/f  collector  noise  formula  in  ultrasmall  BJTs  with  the 
experimental  evidence  shown  in  Sec.  V,  the  first  calculation  of  the  fractal  dimension  of 
the  1/f  noise  process  presented  in  Sec.  VI,  and  new  analytical  quantum  1/f  mobility 
fluctuation  calculations  in  semiconductors.  My  collaborators  have  been  E.  Bernard!,  T. 
Chung,  A.  Fdrst-Chung,  L.M.N.  Sastri,  X.  Hu,  Jian  Xu,  and  M.  Leong.  These  results  will 
be  briefly  presented  below. 


II.  GENERAL  SUFFICIENT  CRITERION  FOR  1/f  NOISE 
IN  CHAOTIC  NONUNEAR  SYSTEMS 

In  spite  of  the  practical  success  of  our  quantum  1/f  theory  in  explaining  electronic  T/f 
noise  in  most  high-tech  devices,  and  in  spite  of  the  conceptual  success  of  our  earlier  classical 
turbulence  approach  to  1/f  noise,  the  question  about  the  ultimate  origin  of  nature’s  omnipresent 
1/f  spectra  remained  unanswered.  During  the  last  three  decades,  we  have  claimed  repeatedly 
that  nonlinearity  is  a  general  cause  of  1/f  noise.  Our  new  result  proves  that  nonlinearity 
always  leads  to  a  1/f  spectrum  if  homogeneity  is  also  present  in  the  equation(s)  of  motion. 
Specifically,  let  the  system  be  described  in  terms  of  the  dimensionless  vector  function  Y(x,t) 
by  the  m^^  order  nonlinear  system  of  differential  equation 

0[t,  X,  Y,  dy/dl  dy/dx^  ...dV/dx^,  d^Y/d\^.  d^Y/dx^^ . an^Y/ax^f^ij  =  q  (  1  ) 
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where  the  vector  function  O  may  be  nonlinear  in  any  of  its  arguments.  If  a  number  9  exists 
such  that  Eq.  (1)  implies 

<D[xet,  xx,  Y,  av/x^at,  av/xaxi  ...aY/xaxp,  b^y/x^w, 

B^y/X'^dx^^ . afTiY/X'Tiaxnf^)  =  0  (2) 

for  any  real  number  X,  the  power  spe<"tral  density  of  any  chaotic  solution  for  the  vector  function 
Y  defined  by  Eq.  (1)  is  proportional  to  1/f. 

In  conclusion,  nonlinearity  +  homogeneity  =  1/f  noise.  The  ultimate  cause  of  the 
ubiquitous  1/f  noise  in  nature  is  the  omnipresence  of  nonlinearities  (no  matter  how  weak)  and 
homogeneity.  The  latter  is  finally  related  to  rotational  invariance  and  to  the  isotropy  of  space. 
All  our  four  specific  theories  of  1/f  chaos  in  nonlinear  systems  are  just  special  cases  to  which 
our  criterion  is  applicable.  They  include  our  magneto-plasma  theory  of  turbulence  in  intrinsic 
symmetric  semiconductors  (1966),  our  similar  theory  for  metals  (1971),  the  quantum  1/f 
theory  (pure  quantum  electrodynamics,  1975),  and  the  theory  of  Musha's  traffic  turbulence 
(1989).  A  fifth  application  v^as  developed  in  March  1991  and  concerns  a  one-dimensional 
crystal,  i.e.,  a  chain  of  atoms  with  slightly  anharmonic  interaction  potentials,  which  is 
presented  next. 


III.  QUANTUM  1/F  FLUCTUATION  SPECTRA  FOR  MODE  ENERGY  AND  PHONON 
NUMBER  IN  A  NONLINEAR  CHAOTIC  CHAIN  OF  ATOMS 

Consider  a  chain  of  atoms  in  the  x  direction,  with  a  lattice  constant  b  and  displacements 
qj  from  the  equilibrium  position.  The  equations  of  motion 

md2qn/dt2  =  A[(qn+i-qn)-(qn-qn-l)]  +  B[(qn  +  i -qn)2-(qn-qn-1 + 

+  C[(qn  +  i-qn)^-(qn-qn-l)^]  (3) 

contain  anharmonic  terms  as  long  as  B  and  C  are  different  from  zero.  With  an=qn+rqn  and 
Pn=qn-qn-i.  neglecting  the  second-order  term,  we  obtain 

md2qn/dt2=  [(qn  +  i -qn)-(qn-qn-1  )][A  +  B(qn+i+qn-l)  +  C(an2+anPn+Pn2)].  (4) 
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Going  over  from  finite  differences  to  a  continuum  description,  we  obtain  a  differential  equation 
for  q(x,t) 

m^2q/^t2  =  b292q/9x2[A  +  2Bq  +  3Cb2(^q/^x)2).  (  5  ) 

Performing  a  fourier  transform  with  respect  to  x, 

m32qk/ot2  =  -Ab^k^qk  -  2b2B|k'2qk'qk-k’dk' 

+  3Cb41dk‘|dk"k*k"(k-k’-k")2qk-qk-qk-k--k"-  (6) 

All  integrals  are  from  minus  infinity  to  infinity.  Substituting  qk=u(k,t)exp[ikbt(A/m)'’^2j^ 
in32u/3t2  T  2ikb(A/m)^/23u/3j  _  -2b2B Jk'^uk’Uk-k’dk' 

+  3Cb‘*|dk'Jdk"k'k'‘(k-k'-k”)2uk'Uk“Uk-k’-k"-  (  7 ) 

our  1/f  chaos  criterion  requires  both  nonlinearity  and  homogeneity,  as  well  as  the  presence  of 
chaos  or  of  a  quasichaotic  state.  The  nonlinearity  condition  is  satisfied  unless  B=0  and  C=0, 
while  homogeneity  requires  the  existence  of  two  numbers  p  and  0  such  that  replacing  k  by  Xk 
everywhere  except  in  the  integration  differentials,  and  replacing  t  by  X®t  leaves  the  equation 
multiplied  by  a  general  factor  XP,  i.e.,  formally  invariant.  In  our  case  we  note  that  this 
criterion  is  satisfied  with  d=2  and  e=-1  if  we  neglect  the  third-order  term  bv  setting  C=Q.  Qn 
the  other  ha"  d.  both  in  the  general  case  and  in  the  8=0  case  the  criterion  is  not  satisfied,  except 
for  some  low-freauencv  limiting  case  in  which  all  k  values  and  frequencies  are  so  small,  that 
we  can  neglect  the  term  with  B  and  one  of  the  left  hand  side  terms. 

To  see  how  the  criterion  works  for  C=0.  and  to  verify  that  we  get  indeed  a  1/f  spectrum 
in  this  case,  as  predicted  by  our  criterion,  we  set  X=1/|k|  and  call  |k|t=z 

m32u/3z2  +  2ib(A/m)^  ^23u/32  =  -2b2B  J(k7k)2uk’Uk-k'dk'  (8) 

Substituting  'j(k,t)=k''* v(k,t)  we  get  for  v 

m32v/3z2  +  2ib(A/m)^ ''23v/3z  =  -2b2B  J(k7k)v(k',z)v(k-k',z)dk7(k-k') 

=  -2b2B|k"v(kk",z)v(k-kk",z)dk7(1-k")  (9) 
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We  note  that  k  has  disappeared  from  the  equation  and  is  present  only  as  a  scale  factor  in  the 
arguments  of  v  on  the  rhs.  Therefore  we  can  expect  the  existence  of  solutions  v(k,z)  of  this 
equation  which  do  not  depend  on  the  first  argument.  Such  solutions  exhibit  "sliding  scale 
invariance",  because  t  and  k  or  t  and  x  provide  a  scale  for  each  other,  with  no  other  scale 
present. 

In  certain  conditions,  instabilities  of  a  solution  of  Eq.  (6)  may  generate  chaos,  or  turbulence. 
In  a  sufficiently  large  system  described  by  the  local  dynamical  equation  (6),  in  which  the 
boundary  conditions  become  immaterial,  homogeneous,  isotropic  turbulence,  (chaos)  can  be 
obtained,  with  a  spectral  density  determined  only  by  Eq.  (6).  In  the  absence  of  instability  and 
chaos,  a  certain  type  of  random  stirring  forces  can  generate  a  quasichaotic  stochastic  state 
which  can  also  be  described  with  our  methods  familiar  from  turbulence  theory.  The  stationary 
autocorrelation  function  A(t)  is  defined  as  an  average  scalar  produci,  the  average  being  over  the 
turbulent  ensemble 

A(t)  =  <u(x,t)u(x,t+T)>  =  j<u|<(t)uk(t+x)>dk  =  jU(k,z)dk.  (10) 


Here  we  have  introduced  the  scalar 

U(k,z)  =  <uk(t)uk(t+x)>  (11) 

of  homogeneous,  isotropic  chaos  (turbulence),  which  depends  only  on  lk|  and  z=|klPT,  because 
there  is  nothing  else  in  Eqs.  (6)  and  (7).  All  integrals  are  from  minus  infinity  to  plus  infinity. 
Isotropy  means  here  the  equivalence  of  the  +x  and  -x  directions.  The  chain  of  integro- 
differential  equations  for  the  correlation  functions  of  any  order  must  obey  the  same  sliding- 
scale  invariance  which  we  have  noticed  in  the  fundamental  dynamical  equation  (7)  above. 
Therefore,  in  isotropic,  homogeneous,  conditions,  uk(t)  can  only  depend  on  k  and  z. 
Furthermore,  the  direct  dependence  on  k  must  reflect  this  sliding-scale  invariance,  and  is 
therefore  of  the  form 

Uk(z)  =  |k(-lv(z).  (12) 

Indeed,  only  this  form  insures  that  uk(z)dk  and  therefore  also  the  corresponding  integrals  and 
multiple  convolutions  in  k  space  have  the  necessary  sliding-scale  invariance. 
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According  to  the  Wiener-Khintchine  theorem,  the  spectral  density  is  the  fourier  transform 
of  A(t), 

Su(0  =  Je27cifxA(t)dT  =  (1/f)  ldk'Jdt‘e2ititV-iv(z)  =  ^f,  (13) 


where  we  have  set  fT=t',  k=tk’,  2=kx=k‘t',  and  the  integral 

C  =  Jdk'jdt’e2^it'k'-iv(z)  =  jdk"jdfe2rtit*k"-1  v(k'*)  (14) 

is  independent  of  f.  We  have  defined  the  vector  k’*=t‘k.  This  confirms  indeed  our  criterion. 

The  1/f  spectrum  obtained  by  us  for  the  amplitude  u(t)  carries  over  also  for  the  squared 
amplitude  u2(t).  indeed,  the  autocorrelation  A'(x)  =  <u2(x,t)u2(x,t+T)>  is  given  by 
2A2(t)+A2(0)  if  we  assume  the  amplitude  u(t)  to  be  well  approximated  by  a  Gaussian  process. 
The  Fourier  transform  of  A^(x)  is  the  autoconvolution  of  C/f,  which  is  C2/f,  if  we  interpret  all 
1/f  spectra  as  the  limit  of  f£“^  spectra  for  arbitrarily  small  e.  Therefore,  Su2(f)  =  2C2/f  + 
A2(0)5(f),  where  5(f)  is  the  delta  function.  The  energy  density  and  the  phonon  number  density 
are  both  proportional  to  u^.  This  proves  that  in  this  case  the  energy  density  and  the  phonon 
number  density  are  both  fluctuating  with  a  1/f  spectral  density  if  they  fluctuate  at  all,  i.e.,  if 
the  system  is  either  chaotic,  or  in  a  quasichaotic  state  caused  by  a  suitable  system  of  random 
stirring  forces.  This  confirms  for  the  one-dimensional  case  the  prediction  of  T.  Musha.  In 
three-dimensional  piezoelectric  crystals  similar  1/f  fluctuations  of  the  phonon  number  are 
predicted  by  the  quantum  1/f  theory  as  we  show  in  Sec.,  and  have  been  observed  experimentally 
in  the  Brillouin  scattering  of  light  by  T.  Musha. 


IV.  QUANTUM  1/F  FLUCTUATIONS  IN  QUARTZ  RESONATORS 

According  to  the  general  quantum  1/f  formula.  r“2Sr(f)=2a/Vf  with  a=e2/ric  =1/137 
and  A=2(Aj/ec)2/37t  is  the  quantum  1/f  effect  in  any  physical  process  rate  r.  Setting 
j=dP/dt=P,  where  P  is  the  vector  of  the  dipole  moment  of  the  quartz  crystal,  we  obtain  for  the 
rate  r  of  phonon  removal  from  the  main  resonator  oscillation  mode  by  scattering  on  a  phonon 
from  any  other  mode  of  the  crystal  the  spectral  density 

Sr(f)  =  r24a(AP)2/37:e2c2,  (15) 
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where  (aP)^  is  the  square  of  the  polarization  rate  change  associated  with  the  removal  of  one  of 
the  N  phonons  present  in  the  main  resonator  mode.  To  calculate  it,  we  write  the  energy  W  of  the 
resonator  and  its  change  in  liie  form 

W  =  (N+1/2)h(0  =  p2/2Vxco2:  AW  =  hco  =PaP/Vxco2.  (16) 

Here  x  is  the  susceptibility  and  V  the  volume  of  the  quartz  crystal.  Solving  the  last 
equation  for  aP,  squaring,  and  multiplying  with  the  first,  we  get 

(aP)2  =  hco^Vx/(2N+1),  (17) 

and 

Sr(f)  =  r24ahw3Vx/37re2c2(2N+1}f  =  r2(2t»2x/37tc3f)  (hw  V/ W) .  (18) 

This  result  is  applicable  to  the  fluctuations  in  the  rate  of  a  single,  well-defined  process. 

The  corresponding  frequency  fluctuations  are  given  by 

to-2Sco(f)  =  (2(o3x/37ic3f){licoV/WQ4)(kT/ha))2,  (19) 

wl'.eic  C  .;.w  quality  fade,  d  the  singd  n  ode  quartz  resonator  considered,  T  the  temperature. 

V.  EXPERIMENTAL  CHECKS  ON  COLLECTOR  QUANTUM  1/f  NOISE  IN  BJTS 

V.1  Introduction' 

1/f  noise  in  bipolar  junction  transistors  (BJTs)  was  treated  by  van  der  Zieini'l^l  who 
applied  a  Hooge-type  approach  similar  to  Kleinpenning's  treatment!'^!  of  pn  junctions,  and  used 
experimental  data  to  determine  the  Hooge  constant  which  was  in  turn  compared  with  the 
quantum  1/f  theory.  However,  since  fhe  BJT  is  a  minority  carrier  device,  it  requires  the 
application  of  the  quantum  1/f  (Handel)  equationt^l'l^l  from  the  beginning,  for  the  correct 
interpretation  of  the  number  of  carriers  in  the  denominator  of  the  Langevin  noise  source. 

In  the  most  elementary  modeit®!  of  a  BJT,  the  collector  current  ic  arises  from  minority 
carriers  injected  from  the  emitter  into  the  base,  which  diffuse  across  the  width  Xb  of  the  base 
and  are  then  all  swept  across  the  reverse-biased  collector  junction  by  the  built-in  field  of  the 


junction.  If  we  neglect  the  usually  small  leakage  current  of  the  collector  junction  and  the  small 
fraction  of  the  carriers  recombining  in  the  base,  we  get  for  a  n+pn  BJT; 

Ic  =  AqDn[noB{exp(qVBE/kT))/XB].  (  1  ) 

where  A  is  the  cross  sectional  area  of  the  base.  q=-e  is  the  charge  of  the  minority  carriers  in 
the  base,  Dn  their  diffusion  coefficient  in  the  base.  nB(0)=noBexp(qVBE/kT)  is  the  electron 
concentration  at  the  limit  of  the  emitter  space  charge  region,  Vbe  is  the  applied  base  -  emitter 
voltage,  and  Xb  is  the  width  of  the  base.  The  expression  in  rectangular  brackets  is  the  electron 
concentration  gradient  calculated  with  the  boundary  condition  of  a  vanishing  electron 
concentration  at  the  limit  of  the  collector  space  charge  region.  We  assurrie  the  base  to  be  much 
narrower  than  the  electron 

diffusion  length  Ln=(Dnx),  XB<<Ln,  but  sufficiently  wide  to  avoid  ballistic  electron  transport 
across  the  base.  Usually  Xb  is  a  fraction  of  a  m'cron. 

Quantum  i/f  fluctuations  of  the  collisional  cross  sections  of  the  electrons  in  the  base 
will  yield  fluctuations  of  the  diffusion  constant,  and  of  the  mobility  (5Dn/Dn=Spy|i); 

5lc  =  Aq(8Dn)[noBexp(qVBE/k  i  )/Xb1.  (  2  ) 

The  corresponding  spectral  density  of  fractional  fluctuations  I'^Sic  is 

lc'^<(5lc)^>f  =  Dn‘2<(5Dn)2>  =  4'2<(5p)2>  =  an/fN.  (  3  ) 

In  the  last  step  our  quantum  1/f  equationt^M^l  was  used,  where  N  is  the  number  of  carriers 
which  define  the  scattered,  or  diffused,  current  leaving  the  base  and  emerging  in  the  collector, 
whiie  On  =  uAn  is  the  effective  quanium  1/?  noice  coefficient,  nr  Hooge  constant.  The  number  of 
electrons  N  is  thus  determined  by  the  effective  lifetime  xc  of  the  electrons,  which  will  be 
slightly  lower  than  the  lifetime  in  the  unbounded  collector  material,  due  to  the  collector  lead 
contact  processes,  and  due  to  lateral  surface  recombination.  Indeed,  we  can  write  N  =  xcic/q. 
Thus  we  finally  obtain  the  spectral  density  of  the  collector  current  fluctuations; 

Sic  =  anlcq/(f^c).  (4) 

in  which  tc  is  the  effective  lifetime  of  the  majority  carriers  in  the  collector.  This  expression 
is  simpler,  but  similar  to  the  expression  derived  earlier,  with  the  important  difference  that 
now  we  have  a  lifetime  of  the  carriers  in  the  denominator,  while  before  it  was  the  usually  much 


smaller  diffusion  time  Td^XB^/Dn  of  the  electrons  in  the  base.  Eq.  (4)  also  implies  that  in 
narrow-base  BJTs  of  various  base-widths  an  will  be  constant,  as  in  other  devices,  rather  than 
ap/xd-  In  the  following  section  we  show  that  this  expression  is  in  good  agreement  with  the 
experimental  data  in  BJTs  with  1/f  collector  noise  spectra. 

V.2  Comparison  of  the  Calculated  Results  with  Experimental  Data 
1.  The  effective  quantum  1/f  noise  coefficient, 
or  Hooge  constant,  an 

We  consider  the  following  scattering  processes^^l  in  the  calculation  of  the  quantum  1/f 
noise  coefficient  an: 

a)  Normal  collision  processes  (Impurity  scattering.  Optical  scattering  and  Acoustical 
phonon  scattering). 

b)  Intervalley  scattering;  there  are  two  types  intervalley  processes,  i.e.,  g-processes 
which  include  Umkiapp,  and  f-processes[''‘^l. 

From  these  points  of  view,  we  obtain  the  current  spectral  density  in  the  form  of  Eq.  (4) 

S|c(f)  =  ct  n  [lcQ'^(^'^c)l  =  cr  A  n  [l^q/fftc) ) 

For  the  case  a),  the  normal  collision  processes: 

an  =aAn  =a4Av2/(37tc2)  =4a(kT/morrc2)  =4.69x  1  O' ’ (5) 

where  we  used  a  =  uoCe2/2h  =  1/137  (the  fine  structure  constant),  c  =  3x10®m/sec,  k  = 

1 .38x1  0'22j/k,  T  =  300K,  and  mg  =  9.1x1  O'^lkg. 

For  the  case  b),  the  intervalley  scattering  +  umkiapp  scattering  1/f  noise  (g- 
processes); 

an  =  aAn  =  a4Av2/(3rtc2)  =  a4(KAk/m)2/37:c2  =5.86x10'^,  (6) 

where  Av  =  Ap/m  =  hAk/m,  |Ak|  =  0.8(27t/a).  and  a  =  5.4A  for  SiliconHS]  jhe  g-proceses 
include  Umkiapp  back  to  the  original  B.Z.,  and  the  conduction  effective  massf^2]  m  _  o.26mo 
Comparing  the  an  =  5.86  x  10'^  with  reference  [14]  Fig.  8  where  the  result  was  from 
the  exact  calculation,  we  find  out  that  they  are  pretty  close.  We  would  like  to  point  out  that  op 
s  aintervalley  which  is  only  an  approximation  due  to  the  high  aimervalley  comparing  with 
CCmpurity  3tld  Oacoustic. 


2.  The  experimental  data  compared  to  the  Hooge  parameters  OHn 


i.  In  n^-p-n  bipolar  transistors: 


SLi(f)/^^  =  ttn/fN. 


This  can  be  written  in  terms  of  the  diffusion  constant  ,  D  =  kTjjyq 


So{f)/D2  =  an/fN 


which  yieldst^^l 


^  ctiin  qlcDn  I  I 

Slc(f)  =  f  InL 


f  Wb^  ^N(We 


where  oh  is  the  Hooge  parameter  for  electrons,  f  is  the  frequency,  Dp  is  the  diffusion  constant 
for  electrons,  wb  is  the  width  of  the  transistor  base  region,  N(0)  is  the  electron  concentration 
for  u.nit  length  at  the  emitter  side  of  the  base,  and  N(wb)  is  the  electron  concentration  for  unit 
length  at  the  collector  side  of  the  base. 

If  we  introduce  the  ratio 


N  (  0 )  ^  ytT^Dn^WB 
N(Wb)  ~  D  n/Ws 


(10) 


where  vp  is  the  saturation  velocity  of  the  electrons  in  the  base  region,  and  the  diffusion  time 

Tdn  =  WB2/2Dn  =  1/27rfT  (11) 

where  ft  is  the  upper  cut  off  frequency  of  the  BJT,  then 


^  /n  “Hn  qlc  .  Vn  +  (Dn/WB)  OHn  .  .  Vn+Dn/Ws 

*'<=<'>■  I  =  I  *"'^'"1  D„/Wb  1 


(12) 


ii.  In  p+-n-p  bipolar  transistors: 


where 


(14) 


P{wb)  ■ 


+  .  and  Tdp  =  WB^/2Dp. 

Up 


and  where  vcp  (=10^  cm/s  in  Si)  is  the  saturation  'elocity  of  holes  in  the  base.  For  details  see[ 
],  and  the  following  comparisons  in  table  I,  where  a  refers  to  the  normal  scattering  quantum  1/f 
calculation  and  b  to  the  more  likely  case  of  g-type  intervalley-umkiapp  scattering.  It  is  more 
'■kely  to  have  this  case  since  it  has  a  much  larger  quantum  1/f  effect,  and  will  mask  the  smaller 
contribution  calculated  for  case  a.  Also,  the  experimentally  noticed  strong  increase  in  1/f  noise 
if  the  transistors  are  c^i  from  single-crystals  so  that  the  current  flows  along  an  [100]  -  like 
direction  confirms  this  and  can  not  be  explained  without  the  quantum  1/f  theory.  The 
experimental  values  in  Table  I  are  much  closer  to  the  values  calculated  for  the  case  b  replacing 
■Tnd  with  the  lifetime  in  the  collector,  xc-  This  corresponds  to  the  inclusion  of  a  corrective 
factor  xnd/xc  in  the  derivation  of  the  collector  quantum  1/f  noise  spectral  density.  However,  a 
proportionality  of  S|c  to  fj  is  noticed  in  many  transistors,  contrary  to  our  suggestion  of 
replacing  tnd  with  tc.  The  present  formula  containing  tc  is  applicable  to  ultrasmall  BJTs  with  a 
very  narrow  base  region  for  which  most  of  the  life  time  of  carriers  diffusing  from  the  emitter 
is  spent  in  the  collector.  This  subject  requires  further  theoretical  and  experimental  study  to 
determine  the  exact  limits  of  applicability  of  the  new  formula. 


1  3 


.  NO.1-MRF90  has  fT=  45OMH2 

:  The  experimental  data  vs  Hooge  parameters  :  (  :  BJTs  2-NEC57867  has  fj  =  8GHz' 
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a:  3.76E-25 
b:  4.69E-22 


VI.  QUANTUM  1/F  NOISE  IS  QUANTUM  CHAOS:  FRACTAL  DIMENSION 
OF  QUANTUM  1/F  NOISE 


Quantum  1/f  noise  in  a  physical  quantity  such  as  a  current  j,  a  cross  section,  a  process 
rate,  or  a  kinetic  coefficient,  such  as  the  mobility,  is  represented  by  an  expression  of  the  form 

5j/<j>  =  a  Ilb(k,X)|2cos(ckt+7kx)  (15) 

k:X 

where  the  sum  is  over  all  electromagnetic  modes  labeled  by  their  wave  vector  k  and  polarization 
X,  with  k=lk|.  Eq.  (15)  is  deterministic  and  the  particle-specific  random  initial  phases  are 
present  in  each  term  as  initial  conditions,  like  the  initial  phases  describing  a  turbulent  fluid 
are  present  in  each  Fourier  component  of  the  perfectly  deterministic  velocity  field.  The 
difference  between  the  mentioned  classical  chaos  (turbulence  )  case  and  the  quantum  1/f  chaos 
we  are  introducing  here  becomes  evident  when  we  recall  that  J  is  in  fact  a  probability  current 
density!  Can  we  calculate  a  spectral  density  of  the  fluctuations  5j  of  this  probability  current 
and  claim  that  it  represents  the  expectation  value  of  the  spectral  density  of  the  observed 
quantum  1/f  fluctuations?  Our  rigorous  derivation  in  second  quantization^?-  18-20]  jglls  us  we 
can,  provided  we  divide  the  single-particle  result  by  N  in  the  case  of  bosons  and  N-1  for 
fermions.  Here  N  is  the  number  of  scattered  particles  used  to  define  what  we  call  scattered 
current  j  or  scattering  cross  section  a.  Encouraged  by  this  result,  we  have  undertaken  an  effort 
to  determine  the  fractal  dimension  of  quantum  1/f  noise  given  by  Eq.  (15)  and  restricted  to  an 
observed  frequency  interval  from  cko=  eo=27tfo  to  a=2kF,  where  eo=27i/T  is  determined  by  the 

duration  T  of  the  1/f  nose  measurement.  Our  objective  is  to  determine  the  fractal  dimension  of 
quantum  1/f  noise  and  to  compare  it  with  measurements  of  the  fractal  dimension  of  physical  1/f 
noise  in  HgCdTe  MWIR  photodiodesl^il.  Calculating  the  spectral  density,  replacing  the 
summation  in  Eq.  (15)  by  an  integral,  and  going  back  to  the  original  fluctuations,  we  obtain  an 
equivalent  representation  of  the  fluctuations(?-i^'20] 

A 

5j/<j>=  2(aA)’'2  |cos(et-t-y£)de/VE  (16) 

where  a=1/137  is  Sommerfeld's  fine  structure  constant  and  aA  =  (2a/37i)(Av/c)2  is  the 
bremsstrahlung  coefficient,  or  infrared  exponent,  of  the  process  which  generates  the  current  j, 
Av  being  the  velocity  change  of  the  charged  particles  in  the  process.  Repeating  the  fluctuations 


periodically  outside  the  observational  interval  T,  we  can  represent  Eq.  (16)  through  an 
equivalent  Fourier  series 


M 

5j/<j>=  2aA  ^cos(neot+Yn)/(neo)^'^  =  x{t)  ('•7) 

n=1 


which  is  similar  to  Eq.  (15),  with  the  terms  of  the  same  frequency  grouped  together,  and  the 
terms  lumped  together  in  harmonics  with  random  phases  Vn.  with  M=A/eo.  To  this  expression 

we  have  applied  the  Takens-Grassberger-Procaccia  analysist22],  by  creating  a  time  series 
xj=x(tj)  with  1<i<N.  We  consider  this  series  as  a  one-dimensional  sample  of  the  quantum  1/f 
process.  We  calculate  the  correlation  function  Ci(*^)(r)  which  is  N'2  times  the  number  of  data 
pairs  (xj.xj)  separated  by  a  distance  less  than  r.  Next  we  form  groups  of  d  consecutive  data 
(xi,xj+i ,...xi+d.i).  Considering  them  as  vectors  in  a  d-dimensional  Euclidean  space,  we  again 
calculate 


N 

Cd(M)(r)  =  N-2  X0[r-llxi-xj||l,  (18) 

i:j  =  i 

where  0(x)  is  the  step  function  (zero  for  negative,  1  for  positive,  and  1/2  for  null  arguments) 
Finally  we  plot  the  curves  logioCd^^V)  as  a  function  of  r,  with  d=0,  1,  2,...  as  a  parameter. 
We  notice  (Fig.  1  for  M=1  and  Fig.  2  for  M=2  and  3  where  for  M=2  the  numbers  d=1,2,3,4, 
have  been  replaced  by  a,b,c,d.)  that  the  slope  of  the  curves  is  just  d,  and  increases  as  d  is 
increased  from  curve  to  curve  up  to  M,  the  number  of  terms  in  Eq.  (17).  We  graphed  the  case 
M=1  (a  single  term  in  Eq.  (17)  in  Fig  1  and  M=2  and  3  in 
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Fig.  2.  After  d=M,  the  slope  does  no  longer  increase.  This  maximal  slope  gives  the  fractal 
dimension  of  the  quantum  1/f  process  with  M  terms.  Our  conclusion  at  this  point  is  that 
quantum  1/f  noise  in  a  frequency  interval  corresponding  to  M  terms  in  Eq.  (18)  is  chaotic  with 
a  fractal  dimension  M.  Had  it  been  stochastic  rather  than  chaotic,  the  slope  of  the  curves  would 
have  kept  increasing  indefinitely. 

Experimentally,  Fote,  Kohn,  Fletcher  and  McDonough  found  a  fractal  dimension  of  dM=10 
for  1/f  noise  measured  in  HgCdTe  MW  infrared  detectors  in  the  interval  from  fo=10'2Hz  to 
F=10  Hz.  Noticing  that  F/fo=103=2’°,  we  expect  M=10,  which  yields  also  d|\/i=10.  This  nice 
agreement  between  theory  and  experiment  indicates  that  the  measured  1/f  noise  is  a  form  of 
(deterministic)  chaos,  whereas  the  other  models  which  competed  with  the  quantum  1/f  theory 
incorrectly  described  it  as  a  stochastic  phenomenon.  Furthermore,  this  agreement  represents 
an  independent  verification  of  the  quantum  1/f  theory,  independent  from  tests  based  on  the 
predicted  magnitude  and  spectral  dependence  of  1/f  noise. 
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